Large-eddy simulation of buoyant airflow in an airborne pathogen transmission scenario

Indoor airflow patterns and the spreading of respiratory air were studied using the large-eddy simulation (LES) computational fluid dynamics (CFD) approach. A large model room with mixing ventilation was investigated. The model setup was motivated by super-spreading of the SARS-CoV-2 virus with a particular focus on a known choir practice setup where one singer infected all the other choir members. The room was heated with radiators at two opposite walls in the cold winter time. The singers produced further heat generating buoyancy in the room. The Reynolds number of the inflow air jets was set to Re=2750, corresponding to an air-changes-per-hour (ACH) value of approximately 3.5. The CFD solver was first validated after which a thorough grid convergence study was performed for the full numerical model room with heat sources. The simulations were then executed over a time of t=20 min to account for slightly more than one air change timescale for three model cases: (1) full setup with heat sources (radiators+singers) in the winter scenario, (2) setup without radiators in a summer scenario, and (3) theoretical setup without buoyancy (uniform temperature). The main findings of the paper are as follows. First, the buoyant flow structures were noted to be significant. This was observed by comparing cases 1/2 with case 3. Second, the dispersion of the respiratory aerosol concentration, modeled as a passive scalar, was noted to be significantly affected by the buoyant flow structures in cases 1–2. In particular, the aerosol cloud was noted to either span the whole room (cases 1–2) or accumulate in the vicinity of the infected singer (case 3). Turbulence was clearly promoted by the interaction of the upward/downward moving warmer/cooler air currents which significantly affected the dispersion of the respiratory aerosols in the room. The study highlights the benefits of high-resolution, unsteady airflow modeling (e.g. LES) for interior design which may consequently also impact predictions on exposure to potentially infectious respiratory aerosols.

Effectiveness of RANS in predicting indoor airborne viral transmission: A critical evaluation against LES

We investigate the dispersal of droplet nuclei inside a canonical room of size 10×10×3.2m3 with a four-way cassette air-conditioning unit placed at the center of the ceiling. We use Reynolds averaged Navier–Stokes (RANS) simulations with three flow rates corresponding to air changes per hour (ACH) values of 2.5, 5, and 10. The room setup as well as the operating conditions are chosen to match those of a recent high-fidelity large eddy simulation (LES) study. We use statistical overloading with a total of one million droplet nuclei being initially distributed randomly with uniform probability within the room. Six nuclei sizes are considered ranging in radius from 0.1 to 10μm (166,667 nuclei per size). The simulations are one-way coupled and employ the Langevin equations to model sub-grid motion. The flow and particle statistics are compared against the reference LES simulations, and we find that the RANS k−ɛ realizable model may be used as a computationally cheaper alternative to LES for predicting pathogen concentration in confined spaces albeit, with potentially increased statistical discrepancy. © 2023 Elsevier Ltd

Large-eddy simulations to define building-specific similarity relationships for natural ventilation flow rates

Natural ventilation can play an important role towards preventing the spread of airborne infections in indoor environments. However, quantifying natural ventilation flow rates is a challenging task due to significant variability in the boundary conditions that drive the flow. In the current study, we propose and validate an efficient strategy for using computational fluid dynamics to assess natural ventilation flow rates under variable conditions, considering the test case of a single-room home in a dense urban slum. The method characterizes the dimensionless ventilation rate as a function of the dimensionless ventilation Richardson number and the wind direction. First, the high-fidelity large-eddy simulation (LES) predictions are validated against full-scale ventilation rate measurements. Next, simulations with identical Richardson numbers, but varying dimensional wind speeds and temperatures, are compared to verify the proposed similarity relationship. Last, the functional form of the similarity relationship is determined based on 32 LES. Validation of the surrogate model against full-scale measurements demonstrates that the proposed strategy can efficiently inform accurate building-specific similarity relationships for natural ventilation flow rates in complex urban environments.

Sensitivity of puff dynamics and airborne droplet nuclei distribution to variations in violent expiration events

We use large eddy simulations to investigate the puff and droplet dynamics from violent expiratory events such as coughs and sneezes in the first few seconds following an ejection. For each of the eleven simulations considered, over 60,000 droplets are ejected and individually tracked using the point-particle Euler–Lagrange approach. We test the sensitivity of the puff and droplet dynamics to various parameters including the ejection volume, momentum, and orientation. We also explore the effect of the mouth shape on the aforementioned dynamics by considering elliptical and circular inlet cross-sections. The results from the simulations compare favorably with a recent theoretical framework put forth by Balachandar et al. (2020) in terms of the puff size and propagation velocity. More importantly however, the theory is able to accurately predict the number and size spectra of the potentially virus-laden droplet nuclei that remain airborne within the puff. We observe that the ejection angle and mouth shape do not significantly affect the puff and droplet dynamics. Additionally, we quantify the carrying capacity of the detached puff portions in terms of the number and size spectra of droplets/droplet nuclei suspended within. © 2022 Elsevier Ltd

Development of a high-speed bioaerosol elimination system for treatment of indoor air.

We developed a high-speed filterless airflow multistage photocatalytic elbow aerosol removal system for the treatment of bioaerosols such as severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Human-generated bioaerosols that diffuse into indoor spaces are 1-10 µm in size, and their selective and rapid treatment can reduce the risk of SARS-CoV-2 infection. A high-speed airflow is necessary to treat large volumes of indoor air over a short period. The proposed system can be used to eliminate viruses in aerosols by forcibly depositing aerosols in a high-speed airflow onto a photocatalyst placed inside the system through inertial force and turbulent diffusion. Because the main component of the deposited bioaerosol is water, it evaporates after colliding with the photocatalyst, and the nonvolatile virus remains on the photocatalytic channel wall. The residual virus on the photocatalytic channel wall is mineralized via photocatalytic oxidation with UVA-LED irradiation in the channel. When this system was operated in a 4.5 m3 aerosol chamber, over 99.8% aerosols in the size range of 1-10 µm were removed within 15 min. The system continued delivering such performance with the continuous introduction of aerosols. Because this system exhibits excellent aerosol removal ability at a flow velocity of 5 m/s or higher, it is more suitable than other reactive air purification systems for treating large-volume spaces.

Development and evaluation of a fluidic facemask for airborne transmission mitigation.

Recently, a fluidic facemask concept was proposed to mitigate the transmission of virus-laden aerosol and droplet infections, such as SARS-CoV-2 (COVID-19). This paper describes an experimental investigation of the first practical fluidic facemask prototype, or "Air-Screen". It employs a small, high-aspect-ratio, crossflow fan mounted on the visor of a filter-covered cap to produce a rectangular air jet, or screen, in front of the wearer's face. The entire assembly weighs less than 200 g. Qualitative flow visualization experiments using a mannequin clearly illustrated the Air-Screen's ability to effectively block airborne droplets (∼100 µm) from the wearer's face. Quantitative experiments to simulate droplets produced during sneezing or a wet cough (∼102 µm) were propelled (via a transmitter) at an average velocity of 50 m/s at 1 m from the mannequin or a target. The Air-Screen blocked 62% of all droplets with a diameter of less than 150 µm. With an Air-Screen active on the transmitter, 99% of all droplets were blocked. When both mannequin and transmitter Air-Screens were active, 99.8% of all droplets were blocked. A mathematical model, based on a weakly-advected jet in a crossflow, was employed to gain greater insight into the experimental results. This investigation highlighted the remarkable blocking effect of the Air-Screen and serves as a basis for a more detailed and comprehensive experimental evaluation.

Effects of upper respiratory tract anatomy and head movement on the buoyant flow and particle dispersion generated in a violent expiratory event.

In the wake of the COVID-19 pandemic, interest in understanding the turbulent dispersion of airborne pathogen-laden particles has significantly increased. The ability of infectious particles to stay afloat and disperse in indoor environments depends on their size, the environmental conditions and the hydrodynamics of the flow generated by the exhalation. In this work we analyze the impact of three different aspects, namely, the buoyancy force, the upper airways geometry and the head rotation during the exhalation on the short-term dispersion. Large-Eddy Simulations have been used to assess the impact of each separate effect on the thermal puff and particle cloud evolution over the first 2 s after the onset of the exhalation. Results obtained during this short-term period suggest that due to the rapid mixing of the turbulent puff, buoyancy forces play a moderate role on the ability of the particles to disperse. Because of the enhanced mixing, buoyancy reduces the range and increases the vertical size of the small particle clouds. In comparison to the fixed frame case, head rotation has been found to notably affect the size and shape of the cloud by enhancing the vertical transport as the exhalation axial direction sweeps vertically during the exhalation. The impact of the upper airway geometry, in comparison to an idealized mouth consisting in a pipe of circular section, has been found to be the largest when it is considered along with the head rotation.

Effects of face shield on an emitter during a cough process: A large-eddy simulation study.

Face shield is a common personal protection equipment for pandemic. In the present work, three-dimensional computational fluid dynamic (CFD) method is used to simulate a cough jet from an emitter who wears a face shield. A realistic manikin model with a simplified mouth cavity is employed. A large eddy simulation with a dynamic structure subgrid scale model is applied to model the turbulence. An Eulerian-Lagrangian approach is adopted to model the two-phase flows, with which the droplets are represented by a cloud of particles. The droplet breakup, evaporation, dispersion, drag force, and wall impingement are considered in this model. An inlet velocity profile that is based on a variable mouth opening area is considered. Special attentions have been put the vortex structure and droplet re-distribution induced by the face shield. It is found that the multiple vortices are formed when the cough jet impinges on the face shield. Some droplets move backward and others move downward after the impinging. It is also found that a small modification of the face shield significantly modifies the flow field and droplet distribution. We conclude that face shield significantly reduces the risk factor in the front of the emitter, meanwhile the risk factor in the back of the emitter increases. When the receiver standing in front of the emitter is shorter than the emitter, the risk is still very high. More attentions should be paid on the design of the face field, clothes cleaning and floor cleaning of the emitters with face shields. Based on the predicted droplet trajectory, a conceptual model for droplet flux is proposed for the scenario with the face shield.

Multiphase large-eddy simulations of human cough jet development and expiratory droplet dispersion

Violent respiratory events play critical roles in the transmission of respiratory diseases, such as coughing and sneezing, between infectious and susceptible individuals. In this work, large-scale multiphase flow large-eddy simulations have been performed to simulate the coughing jet from a human's mouth carrying pathogenic or virus-laden droplets by using a weakly compressible smoothed particle hydrodynamics method. We explicitly model the cough jet ejected from a human mouth in the form of a mixture of two-phase fluids based on the cough velocity profile of the exhalation flow obtained from experimental data and the statistics of the droplets’ sizes. The coupling and interaction between the two expiratory phases and ambient surrounding air are examined based on the interaction between the gas particles and droplet particles. First, the results reveal that the turbulence of the cough jet determines the dispersion of the virus-laden droplets, i.e. whether they fly up evolving into aerosols or fall down to the ground. Second, the droplet particles have significant effects on the evolution of the cough jet turbulence;for example, they increase the complexity and butterfly effect introduced by the turbulence disturbance. Our results show that the prediction of the spreading distance of droplet particles often goes beyond the social distancing rules recommended by the World Health Organization, which reminds us of the risks of exposure if we do not take any protecting protocol.

Sensitivity of airborne transmission of enveloped viruses to seasonal variation in indoor relative humidity

In temperate climates, the peak in infection rates of enveloped viruses during the winter is likely heightened by seasonal variation in relative humidity within indoor spaces. While these seasonal trends are established in influenza and human coronaviruses, the mechanisms driving this seasonality remain poorly understood. Relative humidity impacts the evaporation rate and equilibrium size of airborne particles, which in turn may impact particle removal rates and virion viability. However, the relative importance of these two processes is not known. Here we use the Quadrature-based model of Respiratory Aerosol and Droplets to explore whether the seasonal variation in enveloped viruses is driven by differences in particle removal rates or by differences in virion inactivation rates. Through a large ensemble of simulations, we found that dry indoor conditions typical of winter lead to slower virion inactivation than humid indoor conditions typical of summer;in poorly ventilated spaces, this reduction in inactivation rates increases the airborne concentration of active virions, but this effect was important to virion exposure only when the susceptible person was farther than 2 m downwind of the infectious person. On the other hand, the impact of relative humidity on particle settling velocity did not significantly affect the removal or travel distance of virus-laden particles, suggesting that relative humidity is more likely to affect seasonal transmission via inactivation rates than via particle removal.

Modelling aerosol transport and virus exposure with numerical simulations in relation to SARS-CoV-2 transmission by inhalation indoors.

We provide research findings on the physics of aerosol and droplet dispersion relevant to the hypothesized aerosol transmission of SARS-CoV-2 during the current pandemic. We utilize physics-based modeling at different levels of complexity, along with previous literature on coronaviruses, to investigate the possibility of airborne transmission. The previous literature, our 0D-3D simulations by various physics-based models, and theoretical calculations, indicate that the typical size range of speech and cough originated droplets ( d ⩽ 20 µ m ) allows lingering in the air for O ( 1 h ) so that they could be inhaled. Consistent with the previous literature, numerical evidence on the rapid drying process of even large droplets, up to sizes O ( 100 µ m ) , into droplet nuclei/aerosols is provided. Based on the literature and the public media sources, we provide evidence that the individuals, who have been tested positive on COVID-19, could have been exposed to aerosols/droplet nuclei by inhaling them in significant numbers e.g. O ( 100 ) . By 3D scale-resolving computational fluid dynamics (CFD) simulations, we give various examples on the transport and dilution of aerosols ( d ⩽ 20 µ m ) over distances O ( 10 m ) in generic environments. We study susceptible and infected individuals in generic public places by Monte-Carlo modelling. The developed model takes into account the locally varying aerosol concentration levels which the susceptible accumulate via inhalation. The introduced concept, 'exposure time' to virus containing aerosols is proposed to complement the traditional 'safety distance' thinking. We show that the exposure time to inhale O ( 100 ) aerosols could range from O ( 1 s ) to O ( 1 min ) or even to O ( 1 h ) depending on the situation. The Monte-Carlo simulations, along with the theory, provide clear quantitative insight to the exposure time in different public indoor environments.